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ABSTRACT. Recently, the Navier-Stokes-Voight (NSV) model of viscoelastic incompressible fluid 
has been proposed as a regularization of the 3D Navier-Stokes equations for the purpose of direct 
numerical simulations. In this work we prove that the global attractor of the 3D NSV equations, 
driven by an analytic forcing, consists of analytic functions. A consequence of this result is that the 
spectrum of the solutions of the 3D NSV system, lying on the global attractor, have exponentially 
decaying tail, despite the fact that the equations behave like a damped hyperbolic system, rather 
than the parabolic one. This result provides an additional evidence that the 3D NSV with the small 
regularization parameter enjoys similar statistical properties as the 3D Navier-Stokes equations. 
Finally, we calculate a lower bound for the exponential decaying scale - the scale at which the 
spectrum of the solution start to decay exponentially, and establish a similar bound for the steady 
state solutions of the 3D NSV and 3D Navier-Stokes equations. Our estimate coincides with similar 
available lower bound for the smallest dissipation length scale of solutions of the 3D Navier-Stokes 
equations. 
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1. Introduction 

We consider the Navier-Stokes-Voight (NSV) model of viscoelastic fluid which is governed by 
the system of equations 

u t -vAu-a 2 Au t + (u-V)u + Vp = f, (la) 
Vu = 0, (lb) 
u(x,0) =u in {x), (lc) 

in CI = [0,L] 3 C R 3 , equipped with the periodic boundary conditions. u(x,t) represents the velocity 
field, p is the pressure, v > stands for kinematic viscosity, / is the forcing, and finally a is a 
real positive length scale parameter, for which the ratio ^- characterizes the response time that is 
required for the fluid to respond to the applied force. The system £T|) was first studied by Oskolkov, 
who introduced the NSV equations (see (271, EH) as a model of motion of linear, viscoelastic 
fluid. 

Recently, in 0> the 3D Navier-Stokes-Voight equations were suggested as a regularization 
model for the 3D Navier-Stokes equations, where a is considered a small regularization param- 
eter. First, it was recognized that the inviscid (v = 0) version of the NSV system CD) coincides 
with the inviscid simplified Bardina sub-grid scale model of turbulence. The viscous simplified 
Bardina model was introduced and studied in |[20ll (see also IB, and (H for the original Bardina 
model). In Q the viscous and inviscid simplified Bardina models were shown to be globally 
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well-posed. Moreover, it was also shown that the viscous simplified Bardina model has a finite 
dimensional global attractor, and the energy spectrum was investigated in 0. Viewed from the 
numerical analysis point of view the authors of Q proposed the inviscid simplified Bardina model 
(or equivalently the inviscid NSV equations) as an inviscid regularization (because no additional 
viscosity or hyperviscosity is introduced) of the 3D Euler equations, subject to periodic boundary 
conditions. Motivated by this observation the system (Q]) was also proposed in Q as a regulariza- 
tion, of the 3D Navier-Stokes equations for the purpose of direct numerical simulations for both 
the periodic and the no-slip Dirichlet boundary conditions. 

The addition of the —a 2 Au t term has two main effects. First, it regularizes the equation in a 
way that the three-dimensional system CO becomes now globally well -posed (see Q, |[27l ). On 
the other hand, as was noted in ifTTl . it changes the parabolic character of the original Navier- 
Stokes equations. Therefore, one does not observe any immediate smoothing of the solutions, as 
expected in parabolic PDEs. We also remark that this type of inviscid regularization has been re- 
cently used for the two-dimensional surface quasi-geostropic model |[T8l . In particular, necessary 
and sufficient conditions for the formation of singularity were presented in terms of regularizing 
parameter. 

The long-time dynamics of the system (Q]) has been studied in lfT6l and ifTTll . where the existence 
of the finite-dimensional global attractor of the system has been established. Moreover, upper 
bounds for the number of determining modes, and the fractal dimension of the global attractor of 
the 3D NSV model where derived in ifTTll . In particular, it was shown that the attractor lies in the 
bounded subset of the Sobolev space H (£1), whenever the forcing term / G L 2 (£2). 

In this work we show that the global attractor of the 3D NSV model consists of the real ana- 
lytic functions, whenever the forcing term / is analytic. The idea is to construct an asymptotic 
approximation v(x,t) to the solution u(x,t) of the system (OQ) satisfying 

lim||v(-,/)-«(-,f)lli?(o) = 0, 

and show that v(x,t) lies in certain Gevrey class - a subspace of the real analytic functions. Func- 
tions belonging to Gevrey regularity class are characterized by the exponential decay of the tail 
of their Fourier coefficients. Our method of the proof - splitting of v(x,t) into higher and lower 
Fourier components, has been used before in the context of the weakly damped driven nonlinear 
Schrodinger equation in ll23l and a model of Benard convection in a porous medium in l24l (see 
also QT]). Recently, the authors of Hi followed the same methods to prove the Gevrey regularity 
of the global attractor of the generalized Benjamin-Bona-Mahony equation. 

An important consequence of our result is that the solutions of the 3D NSV system ([]]) lying on 
the global attractor posses a dissipation range, despite the fact that the equations behave like the 
damped hyperbolic system, rather than the parabolic equation. This fact provides an additional 
evidence that (Q]), with the small regularization parameter a, can indeed be used as a model to 
study the statistical properties of turbulent solutions of the 3D Navier-Stokes equations, a subject 
of ongoing research. 

Finally, we obtain bounds for the exponential decaying length scale, which is related to the 
dissipation length scale, of the general solutions of the NSV system lying on the global attractor. 
The obtained estimate is similar to the bounds for the smallest length scale in the turbulent flow 
that was previously calculated for the solutions of the 3D Navier-Stokes equations in Q. In 
addition, using the techniques introduced in ll25l . we estimate the exponential decaying scale of 
the stationary solutions of the 3D NSV and 3D Navier-Stokes equations. Our bounds coincide 
with those obtained in this paper for the general solutions of the NSV system lying on the global 
attractor, and for those of the 3D Navier-Stokes equations reported in Q. 
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2. Preliminaries 

In this paper we will use the following notations, which are standard in the mathematical theory 
of the Navier-Stokes equations (see, e.g., 0, O, Il30l0 . 

Let £2 = [0,L] 3 . We denote by LP(fl), for 1 < p < °o, and H m (Q) - the usual Lebesgue and 
Sobolev spaces of the periodic functions on £2 respectively. Let J be the set of all vector trigono- 
metric polynomials on the periodic domain Q, and denote 

l/ = j(pej :V-cp = 0, and J (p(jc)rfx = oj. 

We set H, and V to be the closures of 1/ in the L 2 {Q.) and H l (Q.) topology respectively. 

We denote by P a : L 2 — > H - the Helmholtz-Leray orthogonal projection operator, and by A = 
— P a A - the Stokes operator subject to the periodic boundary conditions with domain D(A) = 
(H 2 (Q.)) 3 n V. Observe, that in the space-periodic case 

Au = — P a Au = —Au, for all u G D(A). 
The operator A -1 is a positive definite, self-adjoint, compact operator from H into H. We denote 
by < (^z^j = < ^2 < • • • the eigenvalues of A, repeated according to their multiplicities. The 
eigenvalues Xj satisfy, for some dimensionless constant cq > 0, 

— <^<c / /3 , for; = 1,2,3,.... 

In the periodic case this observation is simple (see, e.g., 0), however, in the general case, this is a 
result of the famous Weyl's formula for the case of the Stokes operator due to Metivier (see, e.g., 

ED, EH, (11). 

For any s G R, we can define the Hilbert spaces V s := D(A S / 2 ) with the inner product and norm 

{u,v) s = Uj-vj\j\ 2s , \u\ 2 = (u,u) s , 

for every «,v£ V s , where uj,vj are the corresponding Fourier coefficients of u and v respectively. 
Note, that Vo = H. We will denote the corresponding inner product and norm in H by (-, ■) and | • |, 
respectively. Moreover, we denote V = V\ , and the corresponding inner product and norm will be 
written for u,v G V 

((k,v)) = (m,v)i, ||m|| = Mi. 
For any wi,w>2 € V we define the following bilinear form 

B(wi,w 2 ) =Pc{{w\ • V)w 2 ). 

It can be shown (see, e.g., Q, PUl ) that B can be extended to a continuous map B : V x V — > V, 
where V' = V-\ is a dual space of V. In particular, for u,v,w G V, there exists a constant c > 0, 
depending only on H, such that 

| (b(u,v),w) v , I < cA7 3/4 | M | 1/2 ll«ll 1/2 l|v||||w||, (2) 

where (x,y) v , denotes an action of an element x G V on the element of the dual space y G V. 

Finally, using the above definitions, we write the system (Q]) in the following equivalent func- 
tional form 

u t +vAu + a 2 Au t +B(u,u) = f, (3a) 
w(jc,0) =u in (x). (3b) 
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To show that the solution of the problem ([3]) has an analytic asymptotic (in time) approximation, 
we will use the concept of the Gevrey class regularity. For a given x > 0, and r > 0, we define the 
Gevrey space to be 

G r ,:=D{A r l 2 e^' 2 ) = {ueH : \A r / 2 e' A> ' 2 u\ 2 = £ \uj\ 2 \j\ 2r e 2x ^ < «>}. 

The space is equipped with the corresponding inner product and norm 

( M ,v), T = (A r / 2 e' Al/2 u,A r/2 e zA> ' 2 v) = £ uj- Vj\j\ 2r e 2 ^ , \u\ r , % = \A r l 2 e^'\\, 

for u,v E G{. One can prove that the space of real analytic functions C m (£2) has the following 
characterization 

C m (n) = |J G' x , 

T>0 

for any r > (see, e.g., ED). The concept of the Gevrey class regularity for showing the ana- 
lyticity of the solutions of the Navier-Stokes equations, was first introduced in iTTOl . simplifying 
earlier proofs. Later this technique was extended to the large class of analytic nonlinear parabolic 
equations in flSJ. 

We conclude this section by a few technical propositions that will be used in the proof of 
our main results. First, we will need the following estimates for the nonlinear term. The proof 
of Proposition [T] below, is achieved by standard interpolation estimates using the Gagliardo- 
Nirenberg and Ladyzhenskaya inequalities (see, e.g., 0, QUI ). 

Proposition 1. The bilinear form B(u,u) satisfies: 

(i) If u E V, then B(u,u) E V_i/2> an d 

\B(u,u)\_ 1/2 <cia7 3/4 ||m|| 2 . (4) 

(ii) If u E V3/2. then B(u,u) E H, and 

\B(u,u)\ < C2^j 3 ^ 4 ||m|| I m 1 3/2 ■ (5) 

(iii) For any integer m > 1, ifu E V m+ \, then B(u,u) E V m , and 

Id/ \\ ^ 1 — T/8 1| n 1/4, ,3/4, 1 ,„ 

\B(u,u)\ m < c m Aj ||m|| 7 \uy \u\ m+ i, (6) 
where Ci,C2,c m > are sca/e invariant constants, and c m depends on m. 

Let X > 0, denote by /\ the //-orthogonal projection onto the span of eigenfunctions of A 
corresponding to eigenvalues of the magnitude less then or equal to A. Denote Q\=I — P\. The 
following Poincare-type inequalities hold. 

Proposition 2. Let v E P\G r + l , and v E QxG r +\ Then, 

|v| r +i,T < e zXl "|v|r+i, and |v| r ,T < ^~ 1//2 |v| r +i,T- (7) 

We will also need an estimate for the nonlinear term in the Gevrey space. Similar inequalities 
can be found in |9) (see also Q, iTTOl ). 

Proposition 3. For any X > 0, u, w E G 2 , and v E G\, the following inequality holds 

\(B(u,v),w) l z \ < C{k\ 3/4 |w|j^ 2 |t<|2^ 2 |v|i iT |w|2,T, (8) 
for some scale invariant constant C\ > 0. 
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It is not difficult to prove the following proposition using the Galerkin approximation procedure 
(see HIl). 

Proposition 4. Let seR. Assume that g(t) G L°°([0, T], V s - 2 ), far some < T < oo Then the 
linear problem 

Zt+vAz + a 2 Az t =g(t), z(0)=0, 
has a unique solution z(t) G C([0,r], V s ). In addition, the following estimate holds, 

avv^o 

for all t G [0, T], and d = + a 2 ) 

We will use the following proposition which we state here without a proof. 

Proposition 5. Let q>(?) fte a nonnegative absolutely continuous function on [to, 00 ), far some 
to > 0, satisfying, for all t > to, the inequality 

< -a(p + Z*p 3/2 + c, q>(/ ) = 0. 

Where the positive constant coefficients a,b,c obey the inequality 

'a 

Then for all t > to 

<?(*) <-■ 
a 

Finally, we will need the following Lemma from [15] (see also Q). 

Lemma 1. Leta(t) and b(t) be locally integrable functions on (0,°°) which satisfy for some T > 
the conditions 

1 f'+ T 
liminf — / a(x)dx > 0, 

t^oo t h 

l r'+ T 



bc^<(^ 12 . (10) 



1 r +i 

limsup— / a (%)dl < °°, 

t^oo T Jf 

1 f t+T 
limsup—/ ft + (x)Jx = 0, 

t^oo T J f 



where a = max{— a,0} and b + = max{^,0}. Suppose that §(t) is a nonnegative, absolutely 
continuous function on [0,°°) that satisfies the following inequality, almost everywhere on [0,°°), 

<\>'(t)+a(t)$(t)<b(t). 

Then §(t) — > 0, as t — > oo. 



3. Asymptotic approximation in V m 

The question of global existence and uniqueness of {]} was first studied in ll27l (the inviscid 
case, v = 0, was studied in [31). It was shown that for every u' n G V, the solution of the the system 
(OQ) is globally well -posed, and satisfies u(x,t) G L°°([0,°°), V). In this section we construct an 
asymptotic approximation of the solution of (Q]) in the space V m , for every m > 2. The result can 
be stated as follows. 
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Theorem 1. Let m > 2 be an integer. Consider a solution u(x,t) of the NSV system ([7]), corre- 
sponding to the initial condition u m £ V with the forcing f £ V m _2- Then there exists a function 

vW(*)eL~([0,oo),y m ), 

satisfying 

Iim||«(/)-v (m )(/)|| =0. 

Proof. Let us fix m > 2, and let u m £ V. First, let us write the solution u(t) = v(t) + w(t), where 
v(t) and w(t) satisfy the coupled system 

v,+\Av + a 2 Av, =f-B(u,u), v(0)=0, (11a) 

w t + vAw + a 2 Aw t = 0, w(0) = u'" . (lib) 

This decomposition has been used in ifTTl . First, by using the fact that u(x,t) £ L°°([0,°°), V), 
and applying subsequently the first part of Proposition Q] and Proposition [4] to equation (1 1 lab , we 
conclude that 

v(/)€L~([0,°°),V 3/2 ). (12) 
Next, from equation (lllbl) we immediately get 

\w(t)\ 2 + a 2 \\w(t)\\ 2 < e - vdot (\u in \ 2 + a 2 \\u in \\ 2 ), (13) 

where do = + a 2 ) . Therefore, v(x,t) is an asymptotic (in time) approximation of u(x,t), 
namely 

lirn||w(f)-v(f)|| = lim||w(?)|| =0. 
At the next step, let us consider v^ 2 \x,t) - the solution of the following equation 

v| 2) + vAv (2) +a 2 Av, (2) =f-B(v,v), v (2) (0) = 0. (14) 

According to the second part of Proposition [T] the right-hand side of equation (fl4l is in L°° ( [0, °°) , H) . 
Therefore, applying Proposition |4] we conclude that the unique solution of equation (fl4l satisfies 

v^(t)eL-([0,oo),V 2 ). (15) 

Denote z^ = — v, which satisfies 

+vAz {2) + a 2 Azi 2) =B(u,u-v)+B{u-v,v), z {2) {0)=0. (16) 

According to Proposition |4] equation (fT6l ) has a unique solution z^ 2 \t) £ L°°([0,°°), V3/2). This is 
because « £ L°° ( [0, °°) , V ) , and v satisfies (fl2l . Therefore, we can take an inner product of equation 
( fTBl with z. Using inequality © we get 

^(k (2) (0l 2 + cx 2 |^(0ll 2 )+v|^(0ll 2 = 

= (B(u(t),u(t)-v(t)),z (2) (t)) + (B(u(t)-v(t),v(t)),z^ ) (t)) < 

<cH0-v(0ll(||«(0ll + IK0ll)lk (2) WII< 
<|;ll"(0-v(0ll 2 (ll«(0ll + l|v(0ll) 2 + ^lk (2) (0ll 2 , (iv) 

where the last relation follows from Young's inequality. Finally, we get 

|(k (2) (0l 2 + « 2 ||^ 2 H0ll 2 ) + ^(k (2) (0l 2 +cx 2 ||z( 2 )(0ll 2 )< 

<^K0-v(0ll 2 (ll^)ll + IK0ll) 2 - 
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Using the fact that u(t),v(t) are bounded uniformly in time in the V norm, and that 

IK0-v(f)[| = \\w(t)\\ — >0, as ?^oo, 
we conclude, after applying Lemma Q] that 

lim|| z ( 2 '(0|| =lim||v(0-v (2) (0ll = limllitft) — v( 2 > (0 II = 0. 

We can continue by induction. Fix 2 < n < m, and assume that we have constructed v^\t) G 
L°°([Q,°°),Vj), for j = 2,3, ... ,n — 1, such that for any j 

limllv^'-^CO-vWfOII =lim||«(f)-vW(/)]| =0. (18) 

Let consider the following equation 

v ^ +vAv (") + a 2 Av (n) = /_5( v («-i) )V («-i)) > V W(0)=0. (19) 

Then, according to Proposition |4j and due to the estimates on the nonlinear term of Proposition [T] 
the unique solution v^ n '(t) of the equation (1191 ) satisfies v w (0 GL oo ([0,°°),y n ). Moreover, denote 

z («) =v («) _ v (»-i) ; satisfying 

z (,,) (0)=0. 

Taking the inner product of the last equation with z^ 1 ' (?) and using Proposition Q] and relation (fT8T ) 
we can show by Lemma Q] that 

lim]lz Cn) (0ll =lim||v ( ' | - 2) (f)-v ( ' l " 1) (?)|| =liml|"(0 - v (b) (0II = °> 
t^>°° t^O f->o 

finishing the proof of the Theorem. 

□ 



It can be proved (see ifTTl ) that the solution of the NSV equations (Q~|) satisfies for all t > 

Ht)f< e —^-( K\ 2 + a 2 \\i/ n \\-^ J+^rfe-- ( 20 ) 



Therefore, there exists to, depending on \u m \, \\u m \\ , |/|_i,v,a, and X\, such that for all t > to 

\\u(t)\\ <M i: =^%=. (21) 
ocvv«o 

The following Lemma gives similar bounds for the asymptotic (in time) approximations v^ m >(x,t) 
in the corresponding norms. 

Lemma 2. Let f G V m _2- Consider to > 0, such that the solution of the NSV equations ([7]) satisfies 
the inequality Ml D for all t > fg- the following statements are true: 

(i) The function v(x,t) G L°°([0,o°), V3/2), constructed in Theorem\J} satisfies for all t > to 

\v(t)\ 3 /2<M 3/2 := ^ 7 =(|/|_ 1/2 + c 1 V 3/4 M 2 )- (22) 
OCV y do 

(ii) Thefunction v ly2 \x,t) G L°°([0,°°), V2), constructed in Theorem\J\ satisfies for all t > to 

\v [2) {t)\ 2 <M 2 := -J—d/l+cj^VMart). (23) 

OCV y (jIq 
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(iii) Let m > 2 be an integer. The function v^ m \x,t) G L°°([0,°°), V m ), constructed in Theo- 
rem\Jl satisfies for all t > to 

\v [m \t)\ m <M m := —L—(\f\ m _ 2 + Cm Xf / W l /4 M 3 2 /4 M m ^). (24) 
OCV y ao 

Proof. Recall that v(t) satisfies equation (fTTT >. v^(t) satisfies (fl4l . In general, v^ m \t), for m > 
2, satisfies equation (u"9l . Therefore, the proof of the Lemma is an immediate application of 
Proposition |4j in particular relation ®, and the inequalities of Proposition [T] 

□ 

4. Asymptotic approximation in the Gevrey space G\ 

The results of the previous section show that with a smooth enough forcing the global attractor 
of the system ([]]) lies in C°°(Q), whenever / is C°°(£2). However, our goal is to show that the global 
attractor is real analytic, whenever / is real analytic. For this purpose we use the idea of [23] and 
[24], to construct the asymptotic approximation of the solution of (O in the Gevrey class G\, for 
some x > 0. 

Theorem 2. Let u(x,t) be a solution of the NSV system ([7]), corresponding to the initial condition 
u' n G V with the forcing f G G\ , for some To > 0. Let to > be as in Lemma\2\ then there exists a 
function 

v m (0GL~([? ,°o),G T 2 ), (25) 
for some X > 0, depending only on |/|i x , V, X[ and CH, satisfying 

lim|| M (0-v (B (0|| =0. (26) 

Proof. Let X > to be chosen later. First, consider v^ 2 '(x,f) - an asymptotic approximation of 
u(x,t), which is constructed in Theorem Q] Moreover, according to Lemma [B there exists a 
constant M2 > (see relation (1231). such that 

\v {2) {t)\ 2 <M 2 , V?>? - (27) 

Denote v(t) = i\v^(f), and consider v(t) - a solution of the following equation 

v, + vAv + a 2 Av, + Q x B{v + v, v + v) = /, v(f ) = 0, (28) 

for t > to, where, for notation simplicity, we denoted / = Q\f. The equation (l28l formally looks 
like a projection of the system (Q]) onto the higher wavenumber components, however, the low 
wavenumber modes v of the advection term satisfy a slightly different equation (see also 11261 for 
such a construction for studying data assimilation). Let us denote by 

v (0 (t)=v(t)+v(t), (29) 

for t > to- Our goal is to show first that there exists x > such that v m G G\. Observe, that v is 
just a trigonometric polynomial, and in particular, is analytic. Therefore, we need to show that we 
can choose X large enough, such that v G G\, for some x > 0. Finally, we will show that v (0 (x,t) is 
indeed an asymptotic approximation of u(x,t). 

Note, that in order to prove that the solution of the equation (l28l lies in a Gevrey class of real 
analytic functions we consider the Galerkin procedure to equation (128T ). However, we omit this 
standard procedure, and obtain formal a-priori estimates on the solutions in the relevant Gevrey 
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space norm. Taking formally the inner product of the equation (|28l ) in G\ with v we obtain the 
following inequality 

i|(|vt, + a^|C) + v|v||,< 

< l</,v)i,d + l(B(v,v),v) l>t | + |(B(v,v),v) 1 , t | + |(B{v,w),v),,,| + |(B(v,v),v),,,|. (30) 

Next, we estimate the terms on the right-hand side of 001 ). First, using subsequently the 
Cauchy-Schwartz and Young inequalities, as well as Proposition |2] we get, assuming x < To, 

|(/,v)i,x| < |/V|v| M < ^\f\lx + \\v\l,- (3D 

Next, using Proposition^ Young inequality, and the Poincare-type inequalities of Proposition^ 
we get the following series of estimates for all t > to 

|(5(v,v),v) ljT | < Ci A,] -3 | v| ^ T 2 1 v | ^; T 2 1 v | 2j t < 

^ 5Cj\v\l x \v\ 2 , % y ^ 5C\e^ U1 M]M 2 v 2 

- —^r + - * |vb - m 

|( S (v,v),v) M | < C^U^^Uvh, < CXe ZS/f J2 ^ (33) 

|(5(v,v),v)i, x | <CiXj |v| 2 ; T |v[i, T < ^ 1/4 ^ 3 /4 l v kv ( 34 ) 

|(B(v,v),v)i,,| S <- J^K, £ ^(lC + °X,f ■ ™ 

Let us set x = min{X _1 / 2 ,Xo}. Then, we will choose X large enough satisfying 

/ deM\ /2 M l 2 /2 C,eMx \ V 
max< ' ./ , — - } < -. (36) 

Using the last bounds, we are ready to substitute relations (l3l"l ). (l32l ). (l33l . (l34l) . and (l35l) into 
equation ( f30b - After rearranging the terms, we get 

^(|v|ix + « 2 |v|i T )+^|v|i t < 

Ci /,. |2 , 2| . ,2 \ 3/2 , 5 I/Ii,t , 5CyM\M 2 



Next, using Poincare-type inequality, Proposition |2 and setting d 2 = (1 + a 2 ) _1 , we can write 

^|v| 2 T > ^^\v\l, + a 2 \ v %) > ^(|v| 2 T + a 2 |v| 2 x) . (38) 



Substituting into equation (1371 ) gives us the inequality 

^(lC+^|!,)<-f(l<« + <*'Mi, 



+ (|v|i, + a ; |v|lQ 3/2 + # + 5C;e T 2 - 09) 
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Now we can apply Proposition [5] to the function q>(t) = (|v(f)|j T + a 2 |v(f)|| T ) which is satisfying 
inequality (l39l ). Using relation (fTOb we conclude, that (|v(?) |f T + ot 2 |v(?) || x ) is bounded for all 
t > to, and in particular v(t) G L^Q^o, 00 ),^), whenever the following holds 

C. ( 5\f\j T 5CyM\M 2 \' 2 ' - XV2 

X^Xftf { + 4vX 3/2 ) < 

In order to satisfy the last inequality, we have to choose X large enough, such that 

1 /3 

2 ,1/2,1/2, / C 4 \f\ I, C 5 M\M 2 \ ' 
avX' X> d 2 >^- + —^^ , (40) 

for some absolute constants C^,C=, > 0. For such choice of X we have v(t) G //"([fo, 00 ),^), and 
this proves the first part of the Theorem. 

We are left to show that v a> (x,t) is an asymptotic approximation of the solution u(x,t) of the 
NSV equation (fl. Let z = u - v ffl , and denote z = P\{u - v (2) ), z = Q\u - v. Clearly, by the 
construction and Theorem [H that 

lim\\P x u(t)-v(t)\\ =lim||z(f)|| =0. (41) 

t — >oo f — »oo 

Therefore, to prove (1261 ) we need to show that 

Km\\Q x u{t)-v{t)\\=\im\\z(t)\\=0. 

t — >oo f — >oo 

Observe that £ satisfies the equation 

z,+vAz + a 2 Az t + Q x {B{u,z)+B(z,u)-B{z,z)) =0, z(t ) = Q x u(t ). 
Taking an inner product of the last equation with z we get 

\jt^ '' + a2 ^) +V ^ 2 ■<\(B(z,u),z)\ + \Q k (B(u,z)+B(z,u)-B(z,z),z)\. (42) 

The first summand on the right-hand side of equation (l42l can be estimated as follows. Using Q, 
Proposition |2j and relation (1211) 



for t >tQ. Plugging this inequality into equation (1421) . and using relation (1381 ). we get 

^(kf + « 2 PI| 2 ) + ^(v-^)(|z| 2 + « 2 PI| 2 )<M0, (43) 

where 

HO = l( B («(0^(0) 5 z(0)l + \(. B (z(t),u{t)),z{t))\ + \(B(z{t),z(t)),z(t))\. 
Applying relation (@T|) and using the fact that u is bounded in the V norm, we conclude that 
b{t) — > 0, as t — > oo. Therefore, applying Gronwall's Lemma[T]to equation d43l yields 

lim||£(f)|| =0, 

for X large enough, satisfying 

X>A^) 4 . (44) 



V 

Summarizing, the statement of the Theorem holds for X large enough satisfying relations (l36l) . 
(US) and (EH). 

□ 
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5. Estimating the exponential decaying small scale 

As we have mentioned in the introduction, an additional goal of this research is to provide 
further support for the proposal made in [31 that the NSV system (Q]), with the small regularization 
parameter a, can be used as a numerical model for studying the original Navier-Stokes equations, 
and in particular their statistical properties. Theorem [2] actually states that the global attractor of 
the NSV system consists of real analytic functions u(x,t), whose Fourier spectrum u(k,t) satisfies 
the decay estimate 

\u(k,t)\ <c\k\-' l e~ 

Therefore, following the ideas of Q (see also lfl2l . iTTBl for a different approach), the quantity 
l/X 1 / 1 , can be naturally identified as the exponential decaying length scale, since the exponential 
decay of the spectrum of u is effective only at high wavenumbers satisfying \k\ > A. 1 / 2 . 

In the case of the Navier-Stokes equations the exponential decaying length scale, and similarly 
the radius of analyticity of solutions, can be identified with the smallest effective length scale in 
the turbulent flow (see, e.g., Q, 0, 021, ifTBlO . Classical Kolmogorov theory of turbulence states 
that the smallest effective length scale in the flow is proportional to 
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where 

£ = v(||w|| 2 ), 

is the mean energy dissipation rate, and (•) denotes either the long time average, or the ensemble 
average with respect to the proper invariant probability measure. In Q it was shown that for the 
solution u{t) of the 3D Navier-Stokes equations, as long ||w(f)|| remains bounded uniformly on 
some interval of time [0, T], the smallest length scale of the turbulent flow satisfies 

1~l( 1 -£)\ (45) 



L 

where in the definition of instead of the usual definition of the energy dissipation rate £, the 
authors considered the largest instantaneous energy dissipation rate on the time interval [fi,r], on 
which the solution of the equations remains regular 

e sup = sup v||v(0|| 2 - 

h<t<T 

In the case of the NSV system, similarly to the Navier-Stokes equations, we can define i^sv 
- the exponential decaying length scale. In other words, i^sv is the largest length scale below 
which an exponential decay of the spectrum of the solutions of the NSV system lying on the 
global attractor becomes effective. In this section we would like to derive a lower bound for the 
£nsv> similar to relation (|45T ) for the 3D Navier-Stokes equations. For other estimates on a related 
smallest length scale (via computation of the radius of analyticity of the solutions) of the Navier- 
Stokes equations in 2 and 3 dimensions see lfT2l . |[T3l . and |[T9l (see also O). See also |6l and 
ll9l for other approach to this subject. 

The energy of the NSV system is defined as 

E{t) = | M (f)| 2 + a 2 ||«(0H 2 , 

which satisfies the balance 

~E(t) = -v\\u\\ 2 + (f,u). 
We denote the mean rate of dissipation of energy for the NSV system as 

£ = V(H 2 }, 
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where (•) stands for the long-time average. Moreover, we have the bound 

£<£.»,,, :=vM]. 

In order to find a lower bound for the exponentially decaying length scale of the NSV flow, we 
need to estimate the value of X, from the inequalities d36l ), d4Ql > and (1441 in the proof of Theorem|2j 
since A- -1 / 2 is a lower bound for the radius of analyticity of the solutions of the NSV system lying 
on the attractor, and therefore, 

£nsv>*t 1/2 . 
First, note, that the condition (l44l is satisfied for 

i-^--y^= L ( £ -E)\ (46) 



Moreover, for a small viscosity v and a, we can estimate M%, using the expressions of Lemma |2j 
in the following way 

M 2 ~C 6 



where > is an absolute constant. Therefore, the condition Q6t holds if 

,-1/2 «" 3 „( l *\ l IAT, 



(48) 

M\ Lh sup \LJ 



or, on the other hand if 



Finally, we are left to check when the condition (1401 ) is satisfied. In order to do this, let us assume, 
as it is conventionally done, that / = 0, namely, X is chosen large enough such that the forcing / 
is supported on the modes less than AT 1 . In addition, we assume that X > oc~ 2 , so that d 2 > \oT 2 . 
In that case, the condition (l40l) becomes 



/2 L ll/3 a -2/3M 



and we obtain the estimate 



JT^-JW/ 3 ^) 4 . (49) 



Combining relations (1461) . (1471) . (1481) . and d49l ) we conclude that the exponential decaying length 
scale of the NSV equations satisfies 

£nsv > niin{L,a,L 1 / 3 a 2 / 3 } • (^) 4 . (50) 

Note, that this estimate has the same asymptotic behavior as the estimate of the characteristic 
length scale of the 3D Navier-Stokes equations obtained in Q, without requiring any additional 
assumptions on the regularity of the flow of the system (Q~|). 
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6. Radius of analyticity of stationary solutions 

At the end of the previous section we computed the exponential decaying length scale of the 
NSV model by estimating the radius of analyticity of the functions lying in the global attractor of 
the system. A particular example of the functions lying on the attractor are the stationary solutions 
of the system. The goal of this section is to show that lower bounds for the exponential decaying 
length scale of the stationary solutions of the NSV system are the same as those obtained in the 
last section for the general element of the global attractor. Observe that the NSV equations has 
the same stationary solutions as the 3D Navier-Stokes equations. All calculation in this section 
are formal and can be rigorously justified using the Galerkin approximation procedure. We are 
following the ideas introduced in ll25l . 

The steady state equation of (Q} has the form 

\Au + B{u,u)= f. (51) 
Note the identity, where X is now a dummy variable in the interval [0,c] 

^IMIi.i =2||w||! /2)TJ (52) 

for every time independent function w G G 3 J 2 . Taking an inner product of (f5TT > with A l l 2 e 2xAl12 u, 
we obtain 

V\\u\\l f%i% < \(B(u,u),A 1/2 u)o,x\ + \(f,A 1/2 u)o,x\- (53) 
Let us assume that the forcing / is supported on the first Nf modes. Therefore, we can write 

\ff /i 1/2 \ I ^ 2TAry 2 A ,-l/2| / , \\ ^ 2aA'J /2 )u l/2 ,, ,,2 

\{f,A' u)o^\<e i Ay \{f,u)\<e f Ay V\\u\\ , 

where the last inequality is the result of the energy conservation, and we chose a large a to be 
determined later. Moreover, we can estimate 

c 2 v 

\(B(u,u),A l/2 u)\ Qx < cX\ /4 \\u\\ 2 lx \\u\\ 3/2 ^ < — ||«||* t + -\\u\\j /2z . 



Let us substitute the last two inequalities into equation (1531) to get 

5V|| M ||^ /2 , T <^||,||1 T + , 2 <V 2 V||,|| 2 , (54) 



which we can rewrite 



- C2>l l /2 ||,.l|4 , n„2GNl /2 A A/2n, 



u\\i/ 2 ,z< — j-\\4i,T + 2 e f N/ II "II - (55) 



v 

Applying to the last inequality the identity (1521 we obtain 

^ll«llix<^ll«lltx + 2. 2 <V /2 || M || 2 , 

Denote 

21 3 /2 i3/4 

y(i) = — 1|«|| i,t, F ■=— —vie f Nf \\u\ 
Therefore, y(x) satisfies 

y<y 2 +F 2 < (y + F) 2 . 
Once again, denote, z = y + F, which satisfies 

z<z 2 , 
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therefore 

z(x)=y(?)+F<(z- 1 (0)-'z) ■ 

"Blow-up time" is 

x B > z 1 (0) = —jj- 2 -j± — — > C I > L 

^L|| B ||2 + ^L^«^/^ /2 ||«|| K' \H 2 

where we used the fact that u is a steady state of the system CD, and hence lies in the global 
attractor. As a result, we get u G G\ for all x < x#. Therefore, we showed that the exponential 
decaying length scale of the stationary solutions of the NSV system, satisfies the same bound d50l ) 
as the general solution of the NSV equations lying on the global attractor. Moreover, this bound 
also holds for the smallest length scales of the stationary solutions of the Navier-Stokes equations 
- similar to the general bound obtained in Q. 



7. Conclusions 

We prove that the elements of the global attractor of the 3D NSV equations (Q]) with periodic 
boundary conditions, driven by an analytic forcing, are analytic. A consequence of this result is 
that the solutions of the 3D NSV system lying on the global attractor have exponentially decaying 
spectrum, despite the fact that the addition of the —a 2 Au, term changes the parabolic character of 
the original Navier-Stokes equation, which now starts to behave similar to a damped hyperbolic 
system. 

An important consequence of our result is that the solutions of the 3D NSV system (OQ) lying 
on the global attractor posses a dissipation range - an exponentially decaying spectrum. This fact 
provides an additional evidence that (Q]) with the small regularization parameter a enjoys similar 
statistical properties of the 3D Navier-Stokes equations, a fact that was first suggested in 0. 

Finally, following the ideas of |[25l . we have computed a lower bound of the radius of analyticity 
of the steady state solution of the NSV and Navier-Stokes equations. The bound coincides with 
the one obtained for the general solutions of the system (Q]) lying on the global attractor. 
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